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Abstract 

The Hamiltonian formulation of general relativity on a null surface 
is established in the teleparallel geometry. No particular conditions 
on the tetrads are imposed, such as the time gauge condition. By 
means of a 3+1 decomposition the resulting Hamiltonian arises as a 
completely constrained system. However, it is structurally different 
from the standard Arnowitt-Deser-Misner (ADM) type formulation. 
In this geometrical framework the basic field quantities are tetrads 
that transform under the global SO(3,l) and the torsion tensor. 
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I. Introduction 


The study of asymptotically flat gravitational waves is an important and 
interesting issue in general relativity. It started with the pioneering work 
of Bondi[^, which was subsequently generalized by SachsQ]. It was soon 
realized that the description of gravitational waves on a null surface facili¬ 
tates the characterization of the true, independent degrees of freedom of the 
gravitational held. This characterization may possibly be mandatory to the 
quantization of gravity. Moreover, null surfaces play an important role in the 
study of gravitational radiation. 

Difficulties in working with the dynamics of null surfaces are well known. 
The latter are characterized by the condition = 0. However, if we naively 
impose this condition in Einstein’s equations we spoil the six evolution equa¬ 
tions, since these equations become exempt of second order time derivatives 
and consequently the evolution becomes undetermined. Imposition of the 
above condition in the variation of the Hilbert-Einstein action integral leads 
to nine equations only. Therefore attempts have been made to arrive at a 
well posed characteristic initial value problem. 

The analysis of the initial value problem for asymptotically hat, nonra¬ 
diating space-times is reasonably well understood. Moreover, the Arnowitt- 
Deser-Misner (ADM) Hamiltonian formulation 0 is usually taken as a paradigm 
to the study of the dynamics of spacelike surfaces. In contrast, there does not 
seem to exist a widely accepted formulation of the characteristic initial value 
problem, or of the corresponding Hamiltonian formulation, as we observe 
from the vast literature on the subject. The initial value problem has been 
analysed, for instance, in whereas the Hamiltonian formulation 

has been investigated both in a 2-1-2 decomposition (|^, |^) and in a 3-1-1 
decomposition (Jll], |T^, |3). In particular, the work of refs.|^, is 

developed in the context of Ashtekar variables. While all of these approaches 
add some progress to the understanding of the dynamics of the gravitational 
held on null surafaces, we see that at the present time there does not exist 
a dehnite, irrefutable Hamiltonian formulation which would, accodding to 
Goldberg et. a/.[|T3, display in an isolated form the true degrees of freedom 


and the observables of the theory, in such a way that the dynamics of these 
degrees of freedom is singled out from the dynamics of the remaining held 
quantities. 

In this paper we construct the Hamiltonian for the gravitational held on 
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a null surface in the teleparallel geometry. The analysis of the dynamics 
of spacelike surfaces in this geometry has already been carried out in [|T^]. 
However in that analysis the time gauge condition was imposed in order to 
simplify the considerations. Since we cannot impose at the same time the 
null surface and the time gauge conditions, the problem has to reconsidered 
in a new fashion. 

The analysis of the gravitational held in this geometrical framework has 
proven to be useful, among other reasons because of the appearance of a 
scalar density in the form of a divergence in the Hamiltonian constraint, and 
which is identihed as the gravitational energy density]^. This expression for 
the gravitational energy can be applied to concrete, physical conhgurations of 
the gravitational held (see, for instance, refs. |TB|, 0.00)- In this paper we 
obtain a similar structure. The four constraints of the theory contain each 
one a divergence which altogether constitute a vector density, and which 
strongly suggests a dehnition for the gravitational radiation energy. The 
detailed analysis of this issue is not carried out here. 

One achievement of this long term program is to demonstrate that gen¬ 
eral relativity can be alternatively presented and discussed in the telepar¬ 
allel geometry, without recourse to the Riemann curvature tensor or to the 
Levi-Civita (metric compatible) connection. In this sense, this geometrical 
framework allows an alternative understanding of the gravitational held. 

In section H we present the Lagrangian formulation of the teleparallel 
equivalent of general relativity (TEGR) in a way somewhat diherent from 
what has been presented so far. In the theory is formulated initially 
with a local SO(3,1) symmetry, and in the Hamiltonian analysis, after hxing 
the time gauge condition, it is concluded that in order to arrive at a set of 
hrst class constraints it is necessary to transform the SO(3,l) into a global 
symmetry group. In this paper the symmetry group is taken as the global 
SO(3,l) from the outset. In section HI we present the boundary conditions 
for the tetrad components, assuming that the radiation is due to a localized 
source. In section IV we present in detail the construction of the Hamiltonian, 
obtained by a 3-1-1 decomposition. In the last section we present additional 
comments and point out further developments. 


Notation: spacetime indices ... and SO(3,1) indices a,b,... run from 0 to 
3. In the 3-1-1 decomposition latin indices from the middle of the alphabet 
indicate space indices according to /i = 0,i, a = (0), (i). The tetrad held 
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e“ ^ yields the usual definition of the torsion tensor: T“ = d^e°‘ y — dye°‘ 
The fiat, Minkowski spacetime metric is fixed by rjab = ea/j.eb^g^’’' = (— +++). 


II. The Lagrangian formulation of the TEGR 


In the Lagrangian formulation of the TEGR is presented in terms of 


the tetrad field and a spin connection cu^ab- Both quantities transform under 
the local SO(3,l) group but are not related, not even by the field equations. 
The equivalence of the teleparallel Lagrangian with the Hilbert-Einstein La¬ 
grangian holds provided we require the vanishing of the curvature tensor ten¬ 
sor biJLv{ijj). In the Hamiltonian analysis we conclude that the symmetry 
group must be the global SO(3), and eventually the connection is discarded. 

In this paper we will establish the Lagrangian density in terms of the 
tetrad field only. The symmetry group is the global SO(3,1). The Lagrangian 
density is given by 


L(e) = -keR^^Rbc, (1) 

where k = ^ is Newton’s constant, e = (ief(e“^), Tabc = e.b^ec''Ta^y 

and 


^abc 


_ ^^rjiabc _j_ rjnbac 

4 ^ 




( 2 ) 


Tetrads transform space-time into SO(3,1) indices and vice-versa. The trace 
of the torsion tensor is given by 


T ' rpa 

b — J- ab • 

The tensor is defined such that 

S abcrji _ ^ n-iabcrri , ^ n-iabcrp rjiarji 

abc ^ abc 2 ^ ^ ^ ^ o, ‘ 

The field equations obtained from (1) read 


6L 




0. (3) 
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It can be shown by explicit calculations [^] that these equations yield Ein¬ 
stein’s equations: 


5L 

5e°‘^ 



In order to obtain the canonical formulation we need a hrst order differ¬ 
ential formulation of (1). This is easily obtanained through the introduction 
of an auxiliary held quantity (pabc = —4>acb- The hrst order diherential for¬ 
mulation of the TEGR is described by the following Lagrangian density, 


L(e,0) = keK^^%^abc - 2Tabc) , (4) 

where is dehned in terms of exactly as in terms of 
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4 '" ' " " ^ ' 2 

Variation of the action constructed out of (4) with respect to 0“^^ yields 


(5) 


A — V 

^^abc ^abc i 

which, after some manipulations, can be reduced to 

(6) 

4^abc Tabc • 

The equation above may be split into two equations: 

( 7 ) 

4^a0k d{)€ak : 

(8a) 

4^aik diCd}^ . 

(86) 


Taking into account eq. (7), it can be shown that the second held equation, 
the variation of the action integral with respect to Ca^, leads precisely to (3). 
Therefore (1) and (4) exhibit the same physical content. 

In section IV we will make explicit reference to null surfaces. The theory 
dehned by (1) or (4) describes an arbitrary gravitational held, as there is 
no restriction in the form of a Lagrange multiplier hxing some particular 
geometry. Without going into details we just mention that if we impose the 
condition = 0 in (3) the resulting equation will still have second order 
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time derivatives (note that this equation has one SO(3,1) and one space-time 
index). 

Before closing this section let us make a remark. The theory dehned by 
(1) and (2) has been considered in the literature, in a different context, as 
the translational gauge formulation of Einstein’s general relativity[^]. It is 


argued in this approach that (1) is invariant under local SO(3,l) transforma¬ 
tions up to a total divergence. This divergence is then discarded, from what 
is concluded that (1) exhibits local gauge symmetry. We do not endorse this 
point of view. A careful analysis of this divergence (the last term of eq. (12) 
of [^) shows that in general it does not vanish for arbitrary elements of the 
SO(3,1) group when integrated over the whole spacelike surface. Problems 
arise if the SO(3,l) group elements fall off as const. -|- O(-) when r —>■ oo. 
Therefore the action is not, in general, invariant under such transformations. 
Surface terms play a very important role in action integrals for the gravita¬ 
tional held, so that one cannot arbitrarily add or remove them. Moreover, 
if (1) were actually invariant under the local SO(3,1) group, then the theory 
would have six additional constraints, which would spoil the counting of de¬ 
grees of freedom of the theory (see eqs. (18), (19) and (31) ahead). 


III. The boundary conditions 

In order to guarantee that the space-time of a localized radiating source 
is asymptotically hat we adopt the conditions laid down by Bondiand 
Sachs^] for the metric tensor. The conditions on the tetrads are simply ob¬ 
tained by constructing the tetrads associated with these radiating helds and 
identifying the asymptotic behaviour when r ^ oo. Of course there is an in- 
hnity of tetrads that yield the same metric tensor. However, we will consider 
a typical conhguration and assume the generality of our considerations. For 
simplicity we will consider in detail Bondi’s metric. 

Bondi’s metric is not an exact solution of Einstein’s equations. In terms 
of radiation coordinates (u, r, 6*, 0), where u is the retarded time and r is a 
luminosity distance, Bondi’s radiating metric is written as 

ds^ = —(^e^^ — U‘^r‘^e^'^^dv? — 2e^^dudr — 2Ur‘^e^'^dudO 
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(9) 


+r^ d9‘^ + e sin^d d(f)‘^'^ . 

The metric above is such that the surfaces for which u = constant are null 
hypersurfaces. Each null radial (light) ray is labelled by particular values 
of u, 9 and 0. At spacelike infinity u takes the standard form u = t — r. 
The four quantities appearing in (9), V, U, f3 and 7 are functions of u, r and 
9. Thus (9) displays axial symmetry. A more general form of this metric 
has been given by SachsQ], who showed that the most general metric tensor 
describing asymptotically flat gravitational waves depends on six functions 
of the coordinates. 

The functions in (9) satisfy the following asymptotic behaviour: 


7 


c 

-h ... 

r 


V 

r 


2M 


dd /I / \ 2 


-4d^|lco«- 


(l+ 8 cot^ 6 * 


+. 


1 / dc \ 1 / 3c \ 

where M = M(n, 9) and d = d{u, 9) are the mass aspect and the dipole 
aspect, respectively, and from the function c(n, 9) we define the news function 

dc{u^0) 

du 

One possible realization of this metric tensor in terms of tetrad helds is 
given by 


/ 

—r U e'^ cos9 coscj) 
—r U e'*' cos9 sincj) 
V rU sin9 


(A) “ 2 sin9 coscj) 
gd sincj) 

{^)~^cos9 


0 

r cos9 coscj) 

r cos9 sincj) 

—r e'^ sin9 


0 \ 

—r 6 “'’' sin9 sincj) 
r e~'^ sin9 coscj 

0 / 

( 10 ) 
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From the expression above we obtain the asymptotic behaviour of the 
tetrad components in cartesian coordinates: 


6(0)0 ~ 1 + 0{-) + ... , 

(llo) 

e(o)fc ~ 1 + O(-) + ... , 

(116) 

6(i)0 ~ 0{-) + ... , 

(11c) 

r 


r /In 

^(i)k ^ ^ik “ 1 “ *•• * 

2 r 

(Iki) 


These expressions establish the boundary conditions for the tetrads. As 
a hnal comment, we remark that if we make M = d = 0 in (9), Bondi’s 
metric reduces to the flat space-time metric in radiation coordinates, and so 
does expression (10) for the tetrads. It can be shown that in this case all 
components of the torsion tensor vanish. 


IV. The 3-|-l decomposition 


There are several fundamental differences between the analysis of this section 
and the approach of Goldberg et. a/.|]T^, |^. In the latter, complex valued 
field variables and an orthonormal set of null vectors adapted to a null surface 
are employed. In contrast, we adopt ordinary, real valued tetrads. Never¬ 
theless, the present analysis is conceptually the same as that developed in 
[0, HI- We conclude, however, that it is unecessary to establish a 3-1-1 
decomposition for the tetrads, as it is normally done. The Hamiltonian for¬ 
mulation arises naturally in terms of the four dimensional tetrad held and 
its inverse, as we will see. 

The Hamiltonian formulation is established from the hrst order differen¬ 
tial Lagrangian density (4). Space and time derivatives appear only in Tabc- 
Expression (4) can be rewritten as 


L(e, 0 ) = + ( 12 ) 
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where the dot indicates time derivative, and 


AaOfc _ A 


abc 


eb°ec‘ 


^aij ^ _ 

Thus the momentum canonically conjugated to Cak is given by 

Yiak ^ 

Expression (12) is then rewritten as 


L = e^k - dkCao - 2ke Taij + ke cP^bc • 


jak 


abc 


(13) 


(14) 


In order to establish the Hamiltonian formulation we need to rewrite the 
expression above in terms of Ca/cH"^ and further nondynamical quantities. 


However this is not a trivial procedure. In |]^ the 3+1 decomposition of the 
theory was possible, to a large extent because of the time gauge condition 
e{i) ^ fc = 0. This condition resulted in a tremendous simplihcation of 

the analysis. It is clear that we cannot impose simultaneously the time gauge 
condition and the null surface condition. Therefore the present analysis will 
be totally different from that of . 

The construction can be formally carried out in two steps. First, we 
substitute the Lagrangian held equation (8b) into (14), so that half of the 
auxiliary helds, (paij, are eliminated from the Lagrangian. Second, we should 
be able to express the remaining auxiliary helds, (paok, in terms of the mo¬ 
menta This is a nontrivial step. 

We need to work out the explicit form of H®*^. It is given by 


Y[ak ^ f.^)gOO(^_gkj^a _^aj^k , n„ak 


Oj 


- e 


Oj 


+ 2 e“ (pP Oj) 


0,) + 0, + o,) - + + e^g 




p + *j - g^^^ *j) + 




(16) 






From now on we impose the null surface condition 


= 0 . 

The imposition of this condition at the end of the Legendre transform would 
render inhnities. Denoting (..) and [..] as the symmetric and anti-symmetric 
parts of held quantities, respectively, we can decompose 11“^ into irreducible 
components: 


+ e%n[*^l+ e%n°\ (16) 

where 

= k 0, + 0, - 0,) 

0 , + g^^^ 0 , - g^^^ o,) - 2^*^ o, + , (17a) 

^^k ^ _g0n^^gkj^^ 

n'*"] = k mj + {g^^g^^ - g’^^g^^W mj] ^ k , (18) 

= kep^. (19) 

The crucial observation of this analysis is that only depends on the 
“velocities” Ifl*^! and depend solely on = T°'ij- Therefore we 

can express only six of the “velocity” helds 0“ qj in terms of the momenta 
In order to hnd out which components of can be inverted we 
decompose the latter identically as 

0“ oi = e“ i/jij + e“ aij + e“° Xj , (20) 

with the following dehnitions: 

1 

'4^ij '^ji ~k ^jOi') ) 
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O', 




— O'ji — {4>iOj <PjOi} ) 


— 4>00j 

Substituting (20) in (17a) we find that depends only on 


n<*‘' = A e{2( 

+A“| . (21) 

Therefore if terms like Oij and Xj appear in L, other than in we would 

have difficulties in performing the Legendre transform, because they cannot 
be transformed into any momenta and 11°^ do not depend on them). 

Fortunately, they do not appear. Let us rewrite L given by (14) in terms of 
(15) and (20), assuming from now on that (j)°' ij = T°‘ ij. It is given by 

L = Yi^’^eak + eaodkli^’^-du{eaoYi^^) 

-i-k p (T ■ • — T'^ ■ ■ ■T'^ A 

^ ^ mn-^atj ^ mn^ 13 ' y ^ nk) 

-l'(>.odn“'=+fc|9V‘T“y-e“(9«T‘.,-9^'=T“.,)-2(e“V'-e“V‘)r'd} . 

( 22 ) 

The field (paok appears only in the last line of the expression above. The 
terms that appear together with 11“^ in this line exactly subtract the last 
line of (15). It is possible to check by explicit calculations that the last line 
of (22) can be written as 

. (23) 

We can then proceed and complete the Legendre transform. In the present 
case the latter amounts to expressing ijjik in terms of IfhT. The inversion can 
be made and leads to 
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i-j 1 90m{,90i9jk “ 1 “ 90k9ji)P ^ i^{,90i90k9mn “ 1 “ 90m90n9ik^P r ? ( 2 ^) 


where 


pik ^ ^n(*fc) _ ^ik _ ^ 25 ) 

ke 

Substituting (23) and (24) back in (22) we finally arrive at the primary 
Hamiltonian Hq = pq — L: 

~k ^('’^{dOiQOmQnk ’^dOiQOkQmn) P P ■ (26) 

Since equations (18) and (19) constitute primary constraints, they have 
to be added to Ho, and so the Hamiltonian becomes 


H = Ho + (n°" - ep"^) + 7 ( 7 °° + 4(e,on“'=) . 

(27) 

The quantities aik, Pk and 7 are Lagrange multipliers. 

Next we note that since the momenta {H®'^} are identically vanishing, they 
also constitute primary constraints, which induce the secondary constraints 

= ^ = 0. (28) 

In the proccess of varying H with respect to Cao we only have to consider Hq, 
because variation of the constraints lead to the contraints themselves: 

(n[*^l -kep^^) = (n°^ -kep^)-e^^ (n°* -kepp\ , (29) 

dCao 2 V / 

-—(n°^-A;ep^) = -e“°(n°^-. (30) 

OGaO 
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As we will see, the evaluation of the constraints (7“ according to (28) re¬ 
veals the constraint structure of the Hamiltonian Hq. After a long calculation 
we arrive at 

C- = + ke + \g^^T - g^^T^ 

+l9or{9om9nk - - kc 6^^i^g^^g^^ 

l„0jrpm rpn _|_ njrpO rpm _ r) Okrpj rpn _ r) jkrpQ rpn 1 
\y ^ m-^ raj \ y ran-*- ij ^y ji^ nk ^y ij-^ nk ( 

+ke\e»g„(g^g^ - ]^gmgn,n)P"‘"P“‘ 

+ijo.(jo„,J„s - 59 oi=S„.„)(P”'”7“'' + P“7"“”)} . (31) 

The quantity 7 “*^ appearing in (31) is defined by 

7“‘ = e‘‘'|9«(9“'r,j+j“r‘y)-29“j“TJ,,+ (/'9“+9‘''9“-2j“y'')r«,,} . 

(32) 

We immediately note that (7“ satisfies the relation 

eaoC^ = Ho . (33) 

Therefore we can write the hnal form of the completely constrained Hamil¬ 
tonian: 


H = e„oC“ + a.nnl'‘l - fep'‘) + I3^(n'‘>= - kep>=) + 79 “ + 4(e.on“‘) . 

(34) 

Although Cao appears as a Lagrange multiplier, it is also contained both 
in Hq and in (7“. However, it is possible to check that e^o is a true Lagrange 
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multiplier. By just making use of the orthogonality relations of the tetrads, 
it is possible to verify that the constraints (7“ satisfy the relation 




6C^ 

Sebo 


0 , 


from what we conclude that variation of H given by (34) with respect to Cao 
yields (7“ plus the constraints on the right hand side of (29) and (30). 


V. Comments 


In the last section we have completed the 3+1 decomposition of the La- 
grangian density (4) on a null surface. In this procedure all tetrad and 
metric components are four dimensional quantities. We have not established 
any decomposition for these helds, basically because it was not needed. Since 
the tetrads do not obey any particular gauge condition, the nondynamical 
component e(o)o cannot be identihed with the usual lapse function. 

The hnal form of the Hamiltonian, eq.(34), is written as a sum of the 
constraints of the theory. One major difference between this Hamiltonian 
formulation and the ADM-type formulation is that in the latter the usual 
vector constraint Hi is linear in the momenta^], whereas here both (7*-*^ and 
(7^°^ are linear and quadratic in in general. 

The next step is the determination of the constraint algebra. The alge¬ 
bra of the ten constraints, equations (18), (19) and (31), is expected to be 
quite intricate. The analysis of []^ showed the existence of second class 
constraints. It is likely the same complication arises here. This issue will be 
investigated in the near future. 

As we mentioned in the introduction, one major motivation for the present 
analysis is the establishment of an expression for the gravitational energy- 
momentum vector density. In the present case this expression is restricted 
to conhgurations of the gravitational held that describe gravitational waves. 
Our previous experience on this subject lead us to conclude that the covariant 
gravitational energy-momentum is given by 


= - / d^xd^n 
Jv 


ak 


(35) 
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As before[^], the integral form of the constraint = 0 can be interpreted 
as an energy equation of the type H — E = 0. Expression (35) allows us to 
compute the energy-momentum of the gravitational radiation held for an ar¬ 
bitrary volume of the three-dimensional space. This analysis will be carried 
out in the context of the Bondi and Sachs metrics and presented in detail 
elsewhere. 
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